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19.  Abstract  ^  ,  ,  •  >  * 

Attenuation  apddynamic  moduli  measurements  on  intact  Sierra  White  granite  samples 
were  perfonned^ttilizing  thre^  experimental  techniques:  hysteresis  loop,  resonant  bar,  and 

ultrasonic  spectral  ratio.  Extensional  attenuations  (Q^*4)  and  Young's  moduli  were 
measured  at  0.1  Hz  with  the  hysteresis  loop  and  between  10  and  210  kHz  with  the  resonant 
bar  techniques.  Ultrasonic  P  and  S  wave  velocities  and  spectral  ratio  compressional  (Qp*1)  ~ 
and  shear^Qj'1)  attenuations  were  measured  in  the  100-200  kHz  and  0.6  - 1.1  MHz 
frequency  bands,  then  converted  into  comparable  extensional  attenuations  and  Young's 
modulL;  Compressional  velocities  and  attenuations  were  also  measured  at  40-80  kHz. 

Both  room  dry  and  water  saturated  samples  were  investigated  with  all  three  techniques,  and 
ultrasonic  measurements  at  0.6  - 1.1  MHz  were  conducted  as  a  function  of  confining 
pressure  to  100  MPa. 

The  Young's  moduli  for  room  dry  samples  are  approximately  the  same  (-34  to  35  MPa) 
at  all  frequencies,  and  the  hysteresis  loop  and  resonant  bar  Q^'s  are  small  and  comparable  to 
within  experimental  error  (-0.008  to  0.006).  The  attenuation  mechanism  is  related  to  small 

amounts  of  water  adsorbed  onto  the  surfaces  of  microcracks.  In  contrast,  ultrasonic  Q'l‘s 
for  room  dry  samples  are  much  larger  and  increase  quite  rapidly  with  increasing  frequency 
between  the  three  different  ultrasonic  frequency  bands.  Between  100  and  200  kHz  spectral 
ratio  attenuations  (Qg'1  =  0.024)  are  approximately  3  times  lager  than  resonant  bar 

attenuations  (Qg'1  =  0.0053).  A  scattering  attenuation  mechanism  for  ultrasonic  pulses  is 
suggested  by  these  observations. 

Water  saturation  dramatically  increases  the  hysteresis  loop  and  resonant  bar  Q^'s  by 
factors  of  3  and  7.5,  respectively,  but  has  a  minimal  effect  on  Young's  modulus.  In 
contrast,  the  0.6  - 1.1  MHz  ultrasonic  data  reveal  that  Young's  modulus  nearly  doubles 

while  Qg'1  only  increases  from  0.058  ~ r  °75.  The  0.6  - 1.1  MHz  compressional 
attenuation  Qp-1  actually  decreases  upc  ;  .ration,  from  0.065  to  0.046,  while  the 
40-80  kHz  Qp'1  increases  from  0.020  to  0.064.  The  reduction  of  attenuation  at  the  higher 
frequencies  provides  additional  evidence  of  scattering;  scattering  cross-sections  are  reduced 
when  fluid  is  introduced.  The  large  increase  in  saturated  modulus  reflects  the  unrelaxed 
dynamic  response  of  pore  fluid  at  ultrasonic  frequencies.  Attenuation  is  enhanced  at  lower 
frequencies  by  a  dynamic  viscous  fluid  flow  mechanism.  A  relaxation  peak  in  the 
attenuation  spectrum  below  approximately  10  kHz  is  indicated  by  decreasing  resonant  bar 
Qg^'s  as  a  function  of  increasing  frequency  to  90  kHz. 
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Summary 


In  order  to  clarify  the  frequency  dependence  of  attenuation  between  0.1  and  100  Hz,  there  is  a 
useful  role  for  laboratory  investigations  of  frequency,  strain  amplitude,  pressure,  saturation,  and 
lithology  effects.  If  the  laboratory  properties  of  rock  can  match  seismic  observations  of  both  the 
value  of  Q  and  its  frequency  dependence,  particularly  in  the  1  to  20  Hz  band,  then  seismic 
observations  can  be  interpreted  by  physical  mechanisms  quantified  through  a  series  of  additional 
experiments.  If  not,  appropriate  explanations  must  rely  on  other  physical  realities,  such  as 
contrasts  in  lithology,  anisotropy,  the  presence  of  fractures,  or  scattering  due  to  other 
inhomogeneities. 

In  many  ways  the  best  technique  for  attenuation  and  moduli  determinations  is  the  hysteresis 
loop  technique.  The  seismic  frequency  band  can  be  matched  exactly,  sample  sizes  can  be  adjusted 
to  include  inhomogeneities  and  fractures,  and  the  strain  amplitude  dependence  of  attenuation  to 
very  large  strains  may  be  investigated.  The  first  paper  of  this  report  investigates  the  manner  in 
which  hysteresis  loop  stress-strain  data  is  analyzed  for  Q"1  determination.  An  inconsistency  was 
discovered  in  the  attenuation  definition  Q'1  =  8W/2tiW,  where  8W  is  the  dissipated  energy,  or  the 
area  of  the  hysteresis  loop,  and  W  is  the  maximum  strain  energy  introduced  into  the  sample 
during  one  cycle,  or  the  area  under  the  increasing  load  portion  of  the  loop.  The  inconsistency 
often  results  in  Q'1  that  is  a  factor  of  four  times  less  than  the  equivalent  phase  angle  attenuation, 
Q1  =  tan  <)>.  Where  <j>  is  the  phase  angle  between  stress  and  strain.  The  source  of  the 
inconsistency  is  shown  to  be  the  definition  of  maximum  strain  energy,  which  may  be 
straightforwardly  misapplied  to  the  entire  area  under  a  hysteresis  loop  when  the  experiment  is 
entirely  in  compression.  This  results  in  the  Q1  that  is  four  times  less  than  the  phase  angle  Q'1.  In 
order  to  resolve  the  conflict,  the  origin  of  zero  stress  and  strain  must  always  be  located  at  the 
center  of  the  hysteresis  loop  before  the  definition  Q1  =  8W/2tiW  is  applied.  Accordingly,  the 
resulting  Q'1  must  be  associated  with  the  material  condition  at  the  center  of  the  loop. 

The  second  paper  in  this  report  is  an  investigation  of  the  frequency  dependence  of  attenuation 
and  moduli  for  Sierra  White  granite  utilizing  three  laboratory  techniques:  hysteresis  loop, 
resonant  bar,  and  ultrasonic  pulse.  The  hysteresis  loop  and  resonant  bar  results  for  room  dry 
granite  indicate  that  Young's  modulus  and  extensional  attenuations  (Qg"1  -0.008  to  0.006)  are 
nearly  constant  over  a  broad  frequency  range  including  0.1  Hz  and  from  10  kHz  to  at  least 
200  kHz.  The  attenuation  mechanism  for  room  dry  samples  is  closely  related  to  t^isorbed 
moisture  on  the  surfaces  of  cracks.  For  the  water  saturated  samples,  hysteresis  loop  and  resonant 
bar  extensional  attenuations  increased  dramatically,  by  factors  of  four  to  seven  and  one-half.  The 
corresponding  Young's  modulus  was  either  unchanged  for  the  fundamental  resonant  bar 
frequency  or  decreased  slightly  for  the  hysteresis  loop  experiment.  In  addition,  resonant  bar 
attenuations  in  the  range  10-90  kHz  decreased  with  increasing  frequency,  consistent  with  a 
relaxation  peak  in  the  extensional  attenuation  spectra  at  a  frequency  below  10  kHz.  The 
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attenuation  mechanism  at  these  frequencies  is  local  viscous  flow  of  pore  fluid  in  microcracks  and 
between  adjacent  pore  spaces. 

Ultrasonic  results  in  three  frequency  bands  between  40  kHz  and  1.1  MHz  show  that  although 
the  dry  Young's  modulus  agrees  with  the  resonant  bar  and  hysteresis  loop  data,  attenuation  is 
larger  and  increases  rapidly  with  frequency.  When  the  frequencies  of  the  ultrasonic  and  resonant 
bar  techniques  overlap,  ultrasonic  attenuation  is  larger  by  a  factor  of  3  to  4.  A  dominant 
attenuation  mechanism  at  ultrasonic  frequencies  is  probably  scattering.  The  scattering  attenuation 
mechanism  is  also  combined  with  fluid  losses,  however,  since  Qg'1  and  Qs_1  increase  in  the 
saturated  sample  at  0.6  - 1.1  MHz.  On  the  other  hand,  Qp’1  actually  goes  down  upon  saturation, 
an  observation  that  correlates  well  with  scattering  since  the  cross-sectional  area  of  the  scatterers  is 
reduced  by  the  presence  of  pore  fluid.  Confining  pressure  closes  or  reduces  the  size  of 
microcracks,  increasing  velocities  and  reducing  attenuation.  At  high  pressures  the  dry  and 
saturated  velocities  and  attenuations  converge. 
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The  Definition  and  Calculation  of  Q'1  in 
Hysteresis  Loop  Attenuation  Measurements 

Karl  B.  Coyner 
New  England  Research,  Inc. 

76  Olcott  Drive 

White  River  Junction,  VT  05001 


Abstract.  An  inconsistency  in  the  interpretation  of  experimental  stress-strain  data  for 
attenuation  determinations  is  investigated.  Attenuation  (Q-1)  calculated  from  the  integrated 
areas  of  plotted  hysteresis  loops  depend  on  the  definition  of  maximum  strain  energy.  If 
maximum  strain  energy  is  taken  as  the  entire  area  under  the  loading  segment  of  the 
hysteresis  loop,  the  resulting  Q-1  are  approximately  4  times  less  than  the  Q1  from  the  phase 
angle  <J>  measured  between  cycled  stress  and  strain  (Q1  =  tan  <j>).  This  inconsistency  is 
resolved  by  locating  the  origin  of  zero  stress  and  strain  at  the  center  of  the  hysteresis  loop, 
regardless  of  the  absolute  sign  of  stresses  and  strains,  and  calculating  the  maximum  strain 
energy  during  one  cycle  relative  to  that  origin.  Attenuation  calculated  from  the  ratio  of 
hysteresis  loop  areas  must  be  carefully  interpreted  before  application  to  seismic  wave 
propagation. 


Introduction 

From  a  seismological  perspective  the  hysteresis  loop  technique  probably  provides  the 
most  relevant  and  acceptable  laboratory  determinations  of  rock  attenuation  and  moduli. 
This  is  because  seismic  frequencies  between  approximately  0.01  and  400  Hz  and  strain 
amplitudes  between  10  7  and  1C'3  have  been  attained  (Gordon  and  Davis,  1968;  Walsh  et 
al.,  1970;  Mckavanaugh  and  Stacy,  1974;  Liu  and  Peselnick,  1983;  Jackson  et  al.,  1984; 
Stoll,  1985;  Coyner,  1990).  The  technique  involves  direct  measurements  of  the  stress- 
strain  response  of  rock  during  an  applied  periodic  stress.  This  stress  may  be  uniaxial 
compressional  loading,  torsional  shear,  or  bulk  hydrostatic  compression.  This  method 
provides  critical  experimental  data  that  can  be  used  to  formulate  constitutive  properties  of 
rock.  This  data  can  be  used  in  models  of  high  strain,  nonlinear  attenuation  (Minster  and 
Day,  1986),  and  allows  interpretations  of  fundamental  attenuation  mechanisms  (Spencer, 
198 1).  At  present  there  is  only  a  limited  set  of  laboratory  data  available.  Much  of  it  is 
contradictory  and  perhaps  improperly  presented  as  Q'1  data  because  of  an  inconsistency  in 
the  definition  of  maximum  strain  energy  as  applied  to  hysteresis  loop  analysis. 
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In  the  hysteresis  loop  technique,  attenuation  may  be  calculated  from  either  the  phase 
angle  4>  between  cycled  stress  and  strain  (Spencer,  1981;  Jackson  et  al.,  1984)  or  from  the 
hysteresis  loops  generated  by  plotting  the  stress  versus  strain  data  (Gordon  and  Davis, 
1968;  Walsh  et  al.,  1970;  McKavanagh  and  Stacey,  1974;  Brennan,  1981;  Liu  and 
Peselnick,  1983).  In  the  former,  the  phase  angle  is  determined  by  signal  analysis  of  stress- 
strain  cycles  obtained  with  the  aid  of  computer-based  data  collection  and  analysis 
algorithms.  The  attenuation  is  given  by  the  relation  Q1  =  tan  <j>.  In  the  latter,  the  hysteresis 
loops  are  plotted.  The  attenuation  is  defined  as  Q'1  =  SW  /  27tW,  where  8W  is  the 
dissipated  strain  energy  and  W  is  the  maximum  strain  energy  introduced  into  the  sample 
during  one  stress-strain  cycle.  The  area  of  the  hysteresis  loop  is  proportional  to  8W  and 
the  area  under  the  increasing  load  portion  of  the  loop  is  proportional  to  W.  Q*1  is  calculated 
by  integrating  these  respective  areas  of  the  hysteresis  loop.  Note  that  absolute  areas,  or 
strain  energies,  are  not  important  since  a  ratio  is  calculated. 

Although  these  two  data  reduction  techniques  appear  straightforward,  a  thorough 
comparison  has  apparently  not  been  made  in  previous  experimental  investigations. 
Presumably,  they  should  duplicate  each  other  because  of  the  equivalency 
Q1  =  tan  <J)  =  8W  /  27tW.  A  closer  examination  of  the  techniques  and  published  results, 
however,  has  revealed  a  caveat  in  the  analysis  of  the  plotted  hysteresis  loops.  The  potential 
result  is  an  incorrect  Q'1  that  is  four  times  less  than  the  true  value  or  that  indicated  by  the 
tangent  of  the  phase  angle. 

The  source  of  this  potential  error  is  the  definition  of  maximum  strain  energy  and  the 
relative  location  of  the  origin  for  zero  stress  and  strain.  The  accepted  definition,  from 
which  the  equivalency  is  derived,  is  appropriate  for  sinusoidal  signals  centered  about  zero 
stress  and  strain,  i.e.,  a  rock  undergoing  both  compression  and  tension,  for  which  the 
origin  is  at  the  center  of  the  hysteresis  loop.  In  this  case  the  maximum  strain  energy 
calculated  for  the  compressive  portion  of  the  loop  is  approximately  l/4th  of  the  overall  peak 
to  peak  energy  function. 

If,  on  the  other  hand,  the  hysteresis  loop  is  displaced  on  the  stress-strain  plot,  the 
maximum  strain  energy  changes.  This  is  exactly  what  happens  in  many  laboratory 
measurements,  which  are  normally  biased  by  a  pre-stress  or  simply  run  over  a  sequence  of 
continuously  positive  stresses  and  strains.  For  these  measurements  the  maximum  strain 
energy  stored  in  a  cycle  is  the  whole  area  under  the  hysteresis  loop,  which  is,  equivalently, 
taking  the  origin  at  the  point  of  lowest  stress  and  strain.  The  maximum  strain  energy  is 
approximately  a  factor  of  4  greater  than  the  above  definition.  This  results  in  a  lower  Q-1  by 
a  factor  of  4.  An  accounting  of  this  factor  of  4  is  significant  in  comparing  various 
experimental  results  that  have  utilized  either  of  the  two  techniques. 
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Hysteresis  Loop  Analysis 


The  time-varying  stress  a(t)  and  strain  e(t)  functions  that  describe  a  linear  anelastic  material 
deformed  by  a  steady-state  sinusoidal  stress  of  frequency  to  and  amplitude  o0  are  derived  to 
be  as  follows 


o(t)  =  o0  sin(cot) 
e(t)  =  Jo0  sin(tot  -  <(>). 


Where  J  is  a  compliance.  Since  the  material  is  anelastic  the  strain  lags  stress  in  time  by  the 
phase  angle,  0.  The  strain  function  may  be  expanded, 

e(t)  =  Jt  c0  sin(tot)  -  J2  a0  cos  (tot) 


to  yield  "in  phase",  sin(tot),  and  "out-of-phase"  cos(tot)  strain  components,  where 

Jt  =  J  cos<(>  and  J2  =  J  sin<)>.  The  ratio  of  the  magnitudes,  or  the  real  and  imaginary  parts  of 

the  strain  function,  is  tan  <|>  =  J2  /  Jj. 

The  hysteresis  loop  is  a  plot  of  the  stress-strain  relationship  described  by  these 
functions  and  is  an  ellipse  with  the  origin  centered  at  zero  stress  and  strain.  A  hysteresis 
loop  that  was  generated  by  introducing  a  phase  angle  of  0.0628  radians  between  two  sine 
waves  is  shown  in  Figure  1.  The  maximum  and  minimum  stresses  and  strains  are  ±  a0 
and  ±  J  a0,  respectively,  and  the  hysteresis  loop  is  traced  out  in  a  clockwise  direction. 

The  phase  angle  and  the  hysteresis  loop  are  both  equivalent  expressions  of  the  same 
anelastic  process  through  which  energy  is  absorbed.  A  dimensionless  measure  of 
anelasticity  is  the  Q  factor,  which  is  an  analogue  of  the  Q  used  for  characterizing  the 
efficiency  of  voltage  transfer  in  electric  circuits.  The  inverse  of  Q  may  be  called  the 
dissipation  factor  (Q1),  two  expressions  for  which  are 

Q'1  =  tan  (j) 
and 

Q-i=6W/27tW 

The  first  expression  is  the  "loss  tangent"  and  refers  to  the  tangent  of  the  phase  angle 
between  stress  and  strain.  In  the  second  expression  relative  attenuation  is  obtained  from 
the  ratio  of  8W,  the  energy  dissipated  during  one  cycle,  to  the  maximum  strain  energy  W 
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introduced  into  the  sample  during  one  cycle. 

The  two  definitions  are  equal  for  the  linear  anelastic  material  described  by  the  sinusoidal 
stress  and  strain  time  functions  given  above.  This  is  easily  shown  by  considering  the 
hysteresis  loop  in  Figure  1.  If  stress  times  strain  is  defined  as  strain  energy,  then  the  areas 
within  the  hysteresis  loop  plot  contain  all  of  the  necessary  information  for  calculating  the 
relative  attenuation  energy  ratio.  The  maximum  energy  W  supplied  to  the  material  is  the 
shaded  area  of  Figure  1,  corresponding  to  the  deformation  between  points  A  and  B,  or  the 
integral 


W  = 


The  dissipated  energy  8W  is  the  cross-hatched  area  contained  within  the  hysteresis  loop 
and  is  defined  in  the  integral  form 


8W 


l 


B 

a 

upper 


d£  + 


/ 


C 

®!ower 


d£ 


The  integrals  can  be  easily  solved  by  noting  from  Figure  1  that  W  is  approximately  the 
area  of  a  triangle  with  sides  Oo  and  J  c0>  and  5W  is  the  area  of  an  ellipse  with  semi-major 
axis  along  the  slope  1/J  and  semi-minor  axis  along  the  slope  -1/J.  The  lengths  of  the  axes 
can  be  calculated  from  the  intersections  with  the  equation  for  the  ellipse  obtained  from  the 
stress-strain  relation.  The  two  integrals  can  therefore  be  geometrically  computed  for  small 
attenuations  as 


and 


W  -  J  i  c02  /  2 
5W  »  7t  tan0  Jj  c02 


Substituting  for  the  ratio  W/8W,  it  is  found  that 
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Q1  =  8W/27CW  =  tan<{> 


and  the  definitions  are  thereby  equivalent.  For  the  hysteresis  loop  in  Figure  1, 

Q-i  =  0.0629  (Q  =  15.9). 

For  nonlinear  materials  and  large  attenuations  the  hysteresis  loop  becomes  nonelliptical 
and  nonsymmetrical,  with  cusped  ends  (McKavanagh  and  Stacey,  1974).  In  this  instance 
it  is  necessary  to  integrate  the  areas  of  the  hysteresis  loop  directly  and  to  calculate  Q-1  from 
the  energy  ratio.  The  phase  angle  between  stress  and  strain  is  not  singularly  defined  in  this 
case  as  the  strain  sinusoid  is  distorted. 

Although  the  analysis  of  hysteresis  loops  is  well-known  (see,  for  example,  Lorrain  and 
Corson,  1970),  an  essential  point  is  that  the  stress  and  strain  functions  are  AC-signals,  that 
is,  centered  about  zero  with  equal  positive  and  negative  excursions  and  a  hysteresis  loop 
centered  on  the  origin.  Consequendy,  the  maximum  strain  energy  ~  J  a02  /  2  is  attained  in 
only  one-half  of  the  overall  peak-to-peak  stress  and  strain  amplitudes  ( 2a0  and  2J  a0, 
respectively). 

If  the  sinusoidal  stress  and  strain  time  functions  are  DC-biased,  i.e.,  offset  so  that  they 
are  continuously  positive,  the  definitions  of  maximum  strain  energy  W,  energy  dissipation 
5W,  and  the  identity  of  the  hysteresis  loop  origin  can  become  somewhat  confusing.  In 
Figure  2,  the  same  hysteresis  loop  as  in  Figure  1  is  replotted  with  the  origin  at  zero  stress 
and  strain.  The  stress  and  strain  time  functions  for  this  hysteresis  loop  have  been  shifted 
by  +a  and  +J  a0,  although  with  the  same  peak-to-peak  amplitudes  of  2a0  and  2J  ct0.  The 
energy  dissipated  during  a  cycle  is  still  the  area  of  the  ellipse,  as  in  Figure  1,  or  Jt  tan<|>  Jj 
a02.  Maximum  strain  energy,  however,  or  the  shaded  area  under  the  loading  portion  of  the 
hysteresis  loop  in  Figure  2,  is 

W  =  (2  o0)  tfJjCTo)/ 2  =  2J1a02 


a  factor  of  4  times  greater  than  before.  Accordingly,  when  8W  and  W  are  substituted  into 
the  relative  attenuation  definition,  the  result  is 

Q-1  =  8W/2rcW  =  (tan<t>)  /  4 


Since  the  DC-bia~  has  not  changed  the  phase  angle  between  stress  and  strain,  the  "loss- 
tangent"  definition  of  attenuation  remains  O'1  =  tan<j>.  The  immediate  and  perplexing  result 
is  that  the  attenuation  factor  calculated  from  the  hysteresis  loop  areas  is  l/4th  of  that 
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calculated  from  the  tangeit  of  the  phase  angle. 

This  apparent  inconsistency  can  be  resolved  by  examining  the  maximum  strain  energy 
W  relative  to  the  origin  of  the  hysteresis  loop  in  Figure  2.  From  the  external  perspective, 
the  sinusoidal  stress  deforming  the  material  generates  maximum  strain  energy  that  is  indeed 
2Jj  o02.  For  the  material,  however,  the  state  of  stress  and  strain  appropriate  for  the 
hysteresis  loop  and  energy  calculation  is  not  at  the  origin  of  zero  stress  and  strain  in  Figure 
2,  but  at  the  center  of  the  hysteresis  loop,  i.e.,  as  in  Figure  1.  Although  stress  and  strain 
are  continuously  positive,  the  material  at  +o0  stress  and  +J  a0  strain  is  being  cyclically 
deformed  by  stress  and  strain  functions  with  amplitudes  of  ±  a0  and  ±  J  a0.  Once  this 
reference  at  the  center  of  the  hysteresis  loop  is  Figure  2  is  adopted,  the  maximum  strain 
energy  is  reduced  to  what  it  was  for  the  hysteresis  loop  in  Figure  1,  Jj  o02  /  2.  The 
attenuation  from  the  relative  energy  attenuation  ratio  and  loss  tangent  then  agree, 

Q-1  =  8W  /  2nW  =  tan<|>. 

An  important  point  is  that  the  attenuation  must  be  associated  with  the  deformed 
condition  of  the  material  at  the  center  of  the  hysteresis  loop.  Therefore,  true  attenuation  at 
zero  stress  and  strain  requires  both  tension  and  compression  (or  both  signs  of  shear)  in  the 
hysteresis  loop.  If,  on  the  other  hand,  the  hysteresis  loop  starts  at  zero  stress  and  strain,  as 
in  Fig.  2,  the  measured  attenuation  is  associated  with  the  material  at  +G0  and  +J  cP. 


Discussion 

Although  the  clarification  of  hysteresis  loop  analysis  for  attenuation  is  straightforward, 
the  inconsistency  has  propagated  through  much  of  the  available  experimental  stress-strain 
data.  Therefore,  previous  results  are  re-evaluated  with  respect  to  this  factor  of  4.  This 
resolves  at  least  one  outstanding  discrepancy  and  increases  (by  a  factor  of  4)  several 
sources  of  Q1  data  for  rock.  This  is  particularly  important  if  attenuations  are  to  be 
interpreted  and  compared  with  other  results  (field  observations,  resonant  bar,  ultrasonic). 

There  are  two  techniques  to  measure  and  interpret  attenuation  from  low-frequency 
stress-strain  data.  Either  the  phase  angle  is  measured  directly  (Spencer,  1981;  Jackson  et 
al.,  1984)  or  else  the  hysteresis  loop  is  plotted  out  and  the  areas  integrated  (Gordon  and 
Davis,  1968;  Walsh  et  al.,  1970;  McKavanagh  and  Stacey,  1974;  Brennan,  1981;  Liu  and 
Peselnick,  1983;  Coyner,  1987).  Hysteresis  loop  integration  of  the  areas  shown  in  Figure 
2  and  discussed  in  the  analysis  section  leads  to  a  Q1  that  is  less  than  the  loss  tangent  Q  1 
(tan<j>)  by  a  factor  of  4.  Therefore,  in  most  instances  the  O'1  obtained  by  analyzing  the 
areas  of  hysteresis  loops  have  to  be  increased  by  a  factor  of  4.  This  is  particularly  true  if 
they  are  to  be  interpreted  and  compared  with  the  Q1  results  of  other  experimental 
techniques  (field  observations,  resonant  bar,  ultrasonic.) 
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Table  1  is  a  tabulation  of  experimental  Q-1  from  previous  low  frequency,  stress-strain 
experiments  on  nominally  dry  and  vacuum  dry  rocks  that  were  measured  with  either  the 
phase  angle  or  plotted  hysteresis  loop  technique.  The  rock  type,  frequency,  strain 
amplitude,  and  "reported"  Q4  are  collected  form  the  various  references  given  in  Table  1. 
The  "corrected"  values  have  been  calculated  using  the  analysis  outlined  above,  i.e., 
reported  Q4's  multiplied  by  a  factor  of  4.  Spencer  (1981)  and  Jackson  et  al.  (1984) 
measured  the  phase  angle  directly  and,  consequently,  their  results  do  not  require  correction 
(References  5  and  7).  For  all  of  the  others  (Table  1)  either  the  plotted  hysteresis  loops  or  a 
statement  of  procedure  indicates  that  maximum  strain  energy  was  integrated  under  the  entire 
loop,  resulting  in  an  underestimated  Q4  by  a  factor  of  4.  One  exception  is  Walsh  et  al. 
(1970),  who  report  their  attenuation  results  as  relative  attenuation,  8W  /  W,  and  not  as  Q*1. 
Nevertheless,  in  order  to  interpret  the  results  as  Q4  the  relative  attenuations  still  need  to  be 
divided  by  8it. 

The  strain  amplitude  dependence  of  attenuation  changes  quite  dramatically  if  large  strain 
amplitude  Q4  are  increased  by  a  factor  of  4.  Gordon  and  Davis  (1968)  found  that 
hysteresis  loop  attenuation  data  at  large  strain  amplitudes  (>1CH)  was  essentially  equal  to 
resonant  bar  attenuations  at  low  strain  amplitudes  (<10'6),  even  though  resonant  bar  data 
indicated  a  rapid  increase  in  attenuation  at  intermediate  strain  amplitudes.  The  correction 
increases  hysteresis  loop  attenuations  by  a  factor  of  4.  This  results  in  a  monotonically 
increasing  attenuation  with  strain  amplitude,  independent  of  experimental  technique. 

This  is  a  far  more  plausible  result  because  it  is  consistent  with  models  of  the  frictional 
attenuation  mechanism  (Walsh,  1966;  Mavko,  1979;  Stewart  et  al.,  1983).  At  low  strain 
amplitudes,  less  than  10'6  strain,  attenuation  is  linear  or  independent  of  strain  amplitude. 

Q  1  's  are  quite  low  (Table  1,  references  5, 6,  and  7).  Above  that  point  attenuation 
increases  rapidly  as  nonlinear  friction  becomes  dominant  along  crack  and  pore  surfaces. 

Q1  may  become  extremely  large  at  strain  amplitudes  above  lfr4,  on  the  order  of  10  to 
several  hundred  (1000/Q,  Table  1,  references  1, 2,  and  3).  In  this  group  the  Q4  for 
microcracked  granites  fall  in  the  range  from  15  to  90  (1000/Q,  references  1, 2, 3,  and  4). 
The  Q4  of  several  hundred  is  for  a  friable,  weathered,  Cedar  City  diorite  (Walsh  et  al., 
1970). 

A  discrepancy  that  is  resolved  when  hysteresis  loop  Q4  are  increased  by  a  factor  of  4  is 
the  unexpected  similarity  of  Q4  data  of  Liu  and  Peselnick  (1983)  compared  to  that  of 
Spencer  (1981).  Both  conducted  hysteresis  loop  experiments  on  cylinders  of  granite  at  low 
frequencies  and  low  strain  amplitudes.  Liu  and  Peselnick,  however,  found  a  small  Q4 
(<2.0)  for  room  drv  Westerly  granite.  This  is  a  saturation  condition  that  contrasts  with  the 
large  Q1  (>500)  observed  by  Spencer  for  vacuum  drv  Oklahoma  granite.  These  similar 
results  seem  unusual  since  the  work  of  Tittmann  (1977)  and  Clark  et  al.  (1980)  has 
underscored  the  substantial  increase  in  Q4  caused  by  the  presence  of  volatiles,  particularly 
water.  Resonant  bar  measurements  by  Coyner  (1987)  on  room  dry  Sierra  White  granite 
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found  Q"1  to  be  approximately  0.125,  similar  to  that  found  by  Winkler  et  al.  (1979)  for 
Sierra  White  and  Tittmann  (1984)  for  Westerly  granite.  Therefore,  the  Liu  and  Peselnick 
data  point  appears  anomalous  because  the  Q4  is  so  small  for  a  room-dry  rock  at  zero 
confining  pressure.  If  the  hysteresis  loop  data  of  Liu  and  Peselnick  are  re-interpreted,  and 
the  Q4  increased  by  a  factor  of  4,  to  Q4  <0.0089,  small  Q4‘s  or  typical  microcracked 
granite  are  therefore  preserved  for  vacuum  conditions  (Spencer,  1981)  or  samples  under 
confining  pressure  (Jackson  et  al.,  1984). 

The  hysteresis  loop  analysis  indicates  the  attenuation  factors  deteimined  for  a  particular 
hysteresis  loop  have  to  be  associated  with  the  stress-strain  status  of  the  material  at  the 
center  of  the  hysteresis  loop.  This  has  several  implications  for  experimental  procedure. 
First,  in  order  to  characterize  the  strain  amplitude  dependence  of  a  rock  material  the  center 
of  all  hysteresis  loops  must  coincide.  This  would  presumably  be  at  the  center  of  the  largest 
amplitude  hysteresis  loop,  so  this  implies  a  large  initial  state  of  stress.  Secondly,  in  order 
to  measure  large  strain  amplitudes  at  zero  or  near-zero  stress  a  negative  stress  or  tension 
would  be  required.  A  preliminary  investigation  of  attenuation  is  tension  has  been  made  by 
Bulau  and  Tittmann  (1985). 

Although  modem  experimental  technique  and  digital  signal  processing  favors  the  direct 
measurement  of  phase  angle  in  hysteresis  loop  experiments,  the  analysis  of  hysteresis  loop 
areas  is  still  important  for  problems  in  nonlinear  anelasticity.  At  large  strain  amplitudes 
friction  distorts  the  strain  signal  and  hysteresis  loops  undergo  a  transition  from  elliptical  to 
cusped  shapes  (Brennan  and  Stacey,  1977).  The  phase  angle  is  no  longer  unique  and 
hysteresis  loops  must  be  integrated  in  order  to  obtain  accurate  attenuations.  In  addition,  the 
shape  of  cusped  loops  may  reveal  the  dynamics  of  frictional  processes. 


Conclusion 

It  has  been  shown  that  direct  stress-strain  calculations  of  attenuation  can  lead  to 
inconsistent  results  when  comparing  phase  angle  with  hysteresis  loop  measurements. 

Q4’s  derived  from  the  areas  of  hysteresis  loops  must  be  multiplied  by  a  factor  of  4  in  order 
to  be  comparable  to  those  from  the  phase  angle,  if  the  maximum  strain  energy  is  taken  as 
the  entire  area  under  the  loading  portion  of  the  hysteresis  loop.  Previous  experimental  data 
must  be  interpreted  and  compared  in  light  of  this  connection  factor.  In  particular,  the  strain 
amplitude  dependence  of  attenuation  is  greater  than  previously  realized,  and  room  dry 
granite  can  have  Q4’s  as  large  as  several  hundred  at  amplitudes  greater  than  10  microstrain. 
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TABLE  1.  Reported  and  Corrected  Hysteresis  Loop  Q'1  Factors. 


Rock 

(Dry) 

Frequency 

(Hz) 

Strain 

(10‘6) 

Mode 

1000*Q'1 

Reported  Corrected 

Reference 

Granite 

.0005-.05 

102-103 

Ext 

15 

60 

1 

Quartzite 

.0005-. 05 

lO^lO3 

Ext 

5 

20 

1 

Granite 
(Cedar  City) 

<.05 

103-104 

Ext 

125-91* 

500-360 

2 

Granite 

(Westerly) 

<.05 

500 

Ext 

22* 

88 

2 

Granite,  basalt, 
and  sandstone 

.003-.  1 

1-10 

Ext 

10 

40 

3 

Basalt 

.001-.5 

1.7 

S 

1.90 

7.62 

4 

Granite 

.001-.5 

1.1 

S 

3.76 

15.0 

4 

Sandstone 

.001-.5 

1.2 

S 

13.3-8.0 

52.6-32.3 

4 

Sandstone,  granite,  .004-.4 
limestone  (vacuum  dry) 

0.1 

Ext 

<2.0 

<2.0 

5 

Granite 

(Westerly) 

.01-1 

.01-.1 

Ext 

<2.0 

<8.9 

6 

Granite 

(PC  >10  MPa) 

.33-.003 

0.6 

S 

2.5-0.5 

2.5-0.50 

7 

References:  1,  Gordon  and  Davis,  1968;  2,  Walsh  et  al.,  1970;  3,  McKavanagh  and  Stacey, 
1974;  4,  Brennan,  1981;  5,  Spencer,  1981;  6,  Liu  and  Peselnick,  1983; 

7,  Jackson  et  al.,  1984. 

♦converted  from  reported  relative  attenuations. 
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SINE  WAVE:  .01  CYCLE  PHASE  SHIFT 


CTo 


STRAIN 


Figure  1.  Hysteresis  loop  derived  from  two  sine  waves  offset  by  0.01  cycles,  or 
0.0628  radians,  centered  at  zero  stress  and  strain.  The  shaded  area  is  the  maximum 
energy  W  introduced  into  a  unit  volume  of  sample  during  one  cycle,  corresponding 
to  the  path  from  A  to  B  on  the  hysteresis  loop.  The  interior  area  of  the  hysteresis 
loop  5W  is  the  energy  dissipated  per  unit  volume  during  one  cycle.  The  attenuation 
is  Q'1  =  5W/2kW  =  tan<p  =  (.0628). 


STRESS 


STRAIN 


Figure  2.  Same  hysteresis  loop  as  in  Figure  1,  but  centered  at  +o0  stress  and  +Jo0  strain. 
The  shaded  area  is  the  maximum  energy  introduced  per  unit  volume  during  the  hysteresis 
loop  cycle,  but  this  gives  an  incorrect  relative  attenuation  (V4Q'1)  for  the  material 
at  stress  -kt0  and  +Jo0  strain,  for  which  this  hysteresis  loop  is  representative. 
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Abstract.  Attenuation  and  dynamic  moduli  measurements  on  intact  Sierra  White  granite 
samples  were  performed  utilizing  three  experimental  techniques:  hysteresis  loop,  resonant 

bar,  and  ultrasonic  spectral  ratio.  Extensional  attenuations  (Qg*1)  and  Young's  moduli  were 
measured  at  0.1  Hz  with  the  hysteresis  loop  and  between  10  and  210  kHz  with  the  resonant 
bar  techniques.  Ultrasonic  P  and  S  wave  velocities  and  spectral  ratio  compressional  (Qp*1) 
and  shear  (Q^1)  attenuations  were  measured  in  the  100-200  kHz  and  0.6  -1.1  MHz 
frequency  bands,  then  converted  into  comparable  extensional  attenuations  and  Young's 
moduli.  Compressional  velocities  and  attenuations  were  also  measured  at  40-80  kHz. 

Both  room  dry  and  water  saturated  samples  were  investigated  with  all  three  techniques,  and 
ultrasonic  measurements  at  0.6  - 1.1  MHz  were  conducted  as  a  function  of  confining 
pressure  to  100  MPa. 

The  Young's  moduli  for  room  dry  samples  are  approximately  the  same  (~34  to  35  MPa) 
at  all  frequencies,  and  the  hysteresis  loop  and  resonant  bar  Qg'!’s  are  small  and  comparable  to 
within  experimental  error  (-0.008  to  0.006).  The  attenuation  mechanism  is  related  to  small 
amounts  of  water  adsorbed  onto  the  surfaces  of  microcracks.  In  contrast,  ultrasonic  Q^  's 
for  room  dry  samples  are  much  larger  and  increase  quite  rapidly  with  increasing  frequency 
between  the  three  different  ultrasonic  frequency  bands.  Between  100  and  200  kHz  spectral 

ratio  attenuations  (Qg'1  =  0.024)  are  approximately  3  times  lager  than  resonant  bar 
attenuations  (Qg'1  =  0.0053).  A  scattering  attenuation  mechanism  for  ultrasonic  pulses  is 
suggested  by  these  observations. 

Water  saturation  dramatically  increases  the  hysteresis  loop  and  resonant  bar  Qg^’s  by 
factors  of  3  and  7.5,  respectively,  but  has  a  minimal  effect  on  Young's  modulus.  In 
contrast,  the  0.6  - 1.1  MHz  ultrasonic  data  reveal  that  Young's  modulus  nearly  doubles 
while  Qg*1  only  increases  from  0.058  to  0.075.  The  0.6  - 1.1  MHz  compressional 
attenuation  Qp*1  actually  decreases  upon  saturation,  from  0.065  to  0.046,  while  the 
40-80  kHz  Qp*1  increases  from  0.020  to  0.064.  The  reduction  of  attenuation  at  the  higher 
frequencies  provides  additional  evidence  of  scattering;  scattering  cross-sections  are  reduced 
when  fluid  is  introduced.  The  large  increase  in  saturated  modulus  reflects  the  unrelaxed 
dynamic  response  of  pore  fluid  at  ultrasonic  frequencies.  Attenuation  is  enhanced  at  lower 
frequencies  by  a  dynamic  viscous  fluid  flow  mechanism.  A  relaxation  peak  in  the 
attenuation  spectrum  below  approximately  10  kHz  is  indicated  by  decreasing  resonant  bar 
Qc*1's  as  a  function  of  increasing  frequency  to  90  kHz. 
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Introduction 


Laboratory  determinations  of  dynamic  attenuation  and  moduli  in  rock  reply  on  three 
principal  experimental  techniques:  ultrasonic,  resonant  bar,  and  hysteresis  loop.  They 
differ  substantially  in  both  experimental  arrangement  and  frequency  bandwidth.  Ultrasonic 
experiments  are  typically  performed  at  frequencies  around  1MHz,  with  P  and  S-wave 
attenuations  determined  by  spectral  ratio  analysis  of  broadband  pulses  propagated  through 
both  the  rock  and  a  low  loss  standard  (Toksoz  et  al.,  1979).  At  intermediate  frequencies, 
between  500  Hz  and  100  kHz,  the  resonant  bar  technique  is  utilized  to  measure  attenuation 
from  the  decay  or  spectral  width  of  extensional,  flexural,  and  torsional  resonances  in 
cylindrical  rods  of  material  (Tittmann,  1977,  Winkler,  1979).  At  lower  frequencies,  the 
hysteresis  loop  technique  measures  attenuation  and  modulus  on  samples  subjected  to 
cyclical  stress- strain  paths  at  frequencies  between  0.01  and  500  Hz  (Spencer,  1981; 
Peselnick  and  Outerbridge,  1961;  Gordon  and  Davis,  1968;  Liu  and  Peselnick,  1983). 

In  this  study  attenuation  and  moduli  data  collected  on  Sierra  White  granite  with  the  three 
experimental  techniques  are  compared.  Sierra  White  granite  is  a  fine-grained,  muscovite- 
biotite  granite  with  a  bulk  density  of  2.638  g/cc,  porosity  of  0.8%,  and  a  mean  grain 
diameter  of  0.70  mm.  The  modal  analysis  composition  is  44%  microline,  32%  quartz,  6% 
albite,  1 1%  biotite,  and  5%  muscovite.  The  samples  were  obtained  from  a  single  block  at 
the  Raymond  Quarry  in  Raymond,  California.  Macroscopic  inspection  of  samples  plus 
deformation  and  ultrasonic  P  and  S-wave  velocity  measurements  revealed  homogeneous 
characteristics  with  elastic  anisotropy  of  less  than  5%.  Experiments  were  conducted  on 
several  samples  cored  or  cut  from  the  same  original  quarry  block. 


EXPERIMENTAL  TECHNIQUES  AND  PROCEDURES 
Hysteresis  Loop  Experiments 

The  hysteresis  loop  technique  was  used  to  investigate  extensional  attenuation  and 
dynamic  Young's  modulus  at  low  frequency  in  dry  and  saturated  Sierra  White  granite.  In 
this  technique,  the  strain  response  of  a  sample  subjected  to  a  continuous  sinusoidal  stress  is 
measured  The  technique  was  similar  to  that  employed  by  Spencer  (1981)  and  Liu  and 
Peselnick  (1983).  The  phase  difference  between  strain  and  the  driving  stress  is  one 
definition  of  attenuation,  and  may  be  measured  directly  with  Fourier  phase  analysis. 

Hysteresis  loop  measurements  were  performed  on  a  precisely  ground  cylinder  of  Sierra 
White  granite  that  was  0.15  meters  in  diameter  and  0.51  meters  long.  The  sample  was 
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instrumented  with  three  axial  LVDT  displacement  transducers  mounted  on  two 
circumferential  rings  positioned  with  radial  screws  near  each  end  of  the  sample  (Fig.  1) 

The  LVDT  barrels  were  axially  mounted  in  one  rinb  while  the  LVDT  cores  were  attached  to 
the  second  ring  with  extension  rods.  This  sample  assembly  was  located  in  series  with  a 
precision  load  cell  in  a  servo-controlled  loading  frame.  The  sample  was  exposed  to  the 
atmosphere  at  room  temperature  and  the  sample  column  was  pre-loaded  to  a  constant  axial 
stress  of  0.1  MPa.  A  sinusoidal,  0.1  Hz,  axial  stress  was  then  applied  to  one  end  of  the 
sample  with  a  servo-controlled  ram  operating  in  force  feedback  mode.  A  PC-based  D/A 
converter  was  used  to  generate  the  0. 1  Hz  controlling  signal.  A  voltage  divider  on  the  D/A 
output  controlled  the  amplitude  of  the  force,  thereby  also  controlling  the  strain  amplitude. 
Axial  force  from  the  load  cell  and  sample  displacement  form  the  LVDT  transducers  were 
measured  and  recorded  with  a  PC-based  A/D  converter  as  a  function  of  load  and  strain 
amplitude.  The  data  were  recorded  over  many  cycles  in  order  to  average  noise.  The 
sample  was  measured  both  dry  and  water-saturated  while  exposed  to  the  atmosphere  at 
room  temperature.  A  sheet  of  plastic  film  was  used  to  cover  the  water-saturated  sample  in 
order  to  prevent  drying. 

The  axial  load  and  displacement  data  were  converted  to  axial  stress  and  strain  data  by 
simple  division  with  sample  cross-sectional  area  and  the  ring  separation  distance, 
respectively.  The  three  individual  LVDT  signals  were  averaged.  Extensional  attenuations 
were  calculated  from  the  tangent  of  the  phase  angle  <J)  between  stress  and  strain  determined 
by  FFT  analysis  (Qg*1  =  tan  <(>).  Young’s  moduli  were  calculated  from  the  least-squares  fit 
of  a  straight  line  through  the  hysteresis  loop.  Up  to  several  hundred  stress  cycles  were 
collected  in  order  to  obtain  repeatable  phase  angles.  Low  strain  amplitude  (~10'^)  linear 
anelastic  data  were  used  to  compare  with  the  low  strain  amplitude  data  collected  with  the 
other  techniques. 

The  errors  in  the  load  and  strain  measurements  were  on  the  order  of  1  %  or  less.  The 
phase  angle  calculation  incurred  a  relatively  larger  error  based  on  the  ability  to  obtain 
repeatable  phase  angles,  and  the  corresponding  Q*1  factors  are  precise  to  ±15%. 


Resonant  Bar  Experiments 

Benchtop  resonant  bar  experiments  were  conducted  on  a  slender  square  rod  of  Sierra 
White  granite  0.00635  meters  square  on  a  side  and  0.171  meters  in  length.  The  sample 
was  measured  unjacketed  in  the  room  dry  and  water-saturated  condition  while  exposed  to 
the  atmosphere  at  room  temperature.  The  sample  was  driven  in  extensional  resonance 
modes  with  high  frequency  (1.0  MHz)  compressional  piezoelectric  ceramic  crystals 
epoxied  to  both  sample  ends  (Fig.  2).  One  crystal  was  arbitrarily  identified  as  the 
excitation  source  and  the  second  as  the  receiver.  The  sample  was  supported  in  free  air  on 
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the  bench  with  a  fixture  that  supported  the  sample  at  its  midpoint.  A  lock-in  amplifier  was 
used  to  resonate  the  rod  at  the  fundamental  extensional  resonance  (approximately  10.3 
kHz)  and  at  higher  harmonics,  up  to  the  20th  harmonic  for  the  dry  sample  and  up  to  the  8th 
harmonic  for  the  saturated  sample.  Both  source  and  sensing  piezoelectric  crystals  were 
wired  to  the  amplifier,  one  to  the  reference  sine  wave  output  and  the  second  to  the 
differential  amplifier  input  The  differential  amplifier  input  from  the  sensing  crystal  was 
able  to  "lock-in"  to  the  phase-shifted  sine  wave  frequency  being  used  to  resonate  the  rod. 

A  PC-controlled  D/A  converter  was  used  to  externally  control  the  lock-in  amplifier 
reference  sine  wave  output  to  sweep  through  the  entire  frequency  band  of  interest,  and  an 
A/D  input  was  used  to  measure  the  amplified  signal  from  the  sensing  crystal. 

Extensional  attenuation  was  calculated  by  dividing  the  halfpower  width  (A  f)  of  the 
resonant  peak  by  the  center  frequency  f:  Qg'  *  =  A  f/f.  Extensional  velocity  was  calculated 
by  dividing  the  wavelength  of  the  mode  by  the  center  frequency.  Young's  modulus  was 
calculated  from  the  square  of  extensional  velocity  times  the  density.  An  error  of  less  than 
5%  can  be  associated  with  the  attenuation  and  modulus  from  this  resonant  bar  experiment. 


Ultrasonic  Experiments 

Ultrasonic  velocity  and  attenuation  measurements  were  performed  on  three  different 
Sierra  White  granite  samples  with  three  different  frequency  bandwidth  ultrasonic 
transducers.  The  three  samples  included  the  0.5 1  m  long  cylinder  used  in  the  hysteresis 
loop  experiments,  a  smaller  cylinder  0.076  m  in  diameter  by  0.152  m  long,  and  a  cubic 
block  0.203  m  long  on  each  side.  One  set  of  ultrasonic  transducers  was  capable  of 
propagating  either  P  or  S  waves  in  the  100-200  kHz  frequency  band,  and  a  second  set  in 
the  0.6- 1.1  MHz  frequency  band.  A  set  of  P  wave  transducers  operated  at  40-80  kHz. 

The  low  frequency  transducers  (100-200  kHz  and  40-80  kHz)  were  designed  for  benchtop 
velocity  and  attenuation  measurements  at  very  low  loads  while  the  high  frequency 
transducers 

(0.6-l.lMHz)  were  designed  for  either  benchtop  or  high  confining  and  pore  fluid  pressure 
(to  200  MPa)  measurements. 

Various  ultrasonic  velocity  and  attenuation  measurements  were  made  on  the  Sierra 
White  granite  samples  with  the  different  ultrasonic  transducers.  Measurements  on  the  large 
0.5 1  m  cylinder  were  conducted  in  order  to  compare  with  the  hysteresis  loop  results.  A 
sample  of  this  size  was  necessary  for  the  40-80  kHz  transducer  „  and  benchtop 
measurements  were  made  on  both  the  dry  and  saturated  sample.  The  smaller  0.15  m  long 
cylinder  was  measured  dry  and  saturated  with  the  0.6  - 1.1  MHz  P  and  S-wave  transducers 
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in  both  a  benchtop  configuration  and  as  a  function  of  confining  pressure  to  100  MPa. 
Finally,  the  0.203  m  block  was  measured  in  the  dry  saturation  with  the  100-200  kHz 
transducers.  Measurements  were  also  carried  out  on  aluminum  standards  for  the  spectral 
ratio  analysis. 

The  experimental  apparatus  included  the  specially  constructed  ultrasonic  transducers,  a 
pressure  vessel  system  capable  of  applying  hydrostatic  confining  pressure  and  pore  fluid 
pressure,  and  electronic  equipment  (Fig.  3).  Ultrasonic  transducers  were  developed  and 
constructed  to  produce  reproducible  signals  at  the  various  frequencies.  The  higher 
frequency  transducers  incorporated  both  P  and  S  crystals  so  that  either  signal  could  be 
propagated  during  an  experiment.  At  50  kHz  a  satisfactory  S-wave  signal  could  not  be 
produced  and  only  a  P-wave  crystal  was  used. 

The  electronic  equipment  used  in  the  experiments  include  switches,  a  pulse  generator, 
amplifier,  filter,  and  digital  oscilloscope,  and  a  computer  (Fig.  3).  a  source  and  receiver 
pair  of  like  crystals  (P  or  S)  were  selected,  and  the  source  crystal  was  excited  by  a  fast 
risetime  electrical  pulse  generated  with  a  ranametrics  5055PR  pulser-receiver  or  a  higher 
power  pulse.  The  crystal  produced  a  broadband  ultrasonic  pulse  that  propagated  through 
the  adjacent  metal  plate  to  which  the  crystal  was  epoxied,  through  the  rock  sample  along  the 
core  axis,  through  the  metal  plate  of  the  receiving  transducer,  and  into  the  receiver  crystal. 
The  electrical  signal  produced  by  the  receiver  transducer  was  amplified  and  filtered  by  the 
receiving  section  of  the  Panametrics  5055PR,  then  digitized  by  a  LeCroy  9400  digital 
oscilloscope.  The  signals  were  typically  amplified  by  20  to  40  dB  and  high-pass  filtered. 
The  oscilloscope  displayed  the  signal  for  the  travel  time  pick  and  the  digitized  signal  was 
transferred  to  the  computer  for  subsequent  spectral  attenuation  analysis.  With  this  system  it 
was  possible  to  acquire  well-defined,  broadband  P  and  S  wave  ultrasonic  signals, 
examples  of  which  will  be  provided  in  the  results  section. 

An  interactive  software  routine  was  developed  to  calculate  attenuation  using  a  spectral 
ratio  technique  described  by  Toksoz  et  al.  (1979).  In  the  experimental  procedure,  two 
signals  were  recorded  with  the  equipment  settings  exactly  the  same,  one  through  the  rock 
and  the  second  through  the  aluminum  standard  of  the  same  dimensions.  With  the  computer 
the  digitized  signals  were  displayed  and  windowed  in  the  time  domain  over  approximately 
1.5  to  2  full  cycles.  The  window  edges  were  smoothed  with  a  cosine  taper.  The  two 
windowed  signals  were  decomposed  into  amplitude  versus  frequency  data  utilizing 
standard  Fourier  analysis.  The  individual  amplitude  data  and  the  logarithm  of  the  ratio 
between  them  were  displayed  as  a  function  of  frequency.  A  least  squares  fit  of  the  slope 
was  made  over  a  frequency  band  defined  interactively.  The  inverse  quality  factor  Q"1  was 
calculated  by  multiplying  the  slope,  in  units  of  time,  with  sample  velocity,  and  dividing  by 
the  sample  length  times  pi  (tt). 

This  method  of  measuring  and  calculating  attenuation  requires  several  assumptions. 
Attenuation  is  assumed  to  be  constant  over  the  frequency  band  of  the  signal,  geometrical 


20 


factors  such  as  spreading  and  reflections  are  assumed  constant  for  the  two  signals,  and 
attenuation  of  the  aluminum  is  taken  as  very  small  relative  to  the  rock.  An  important 
variable  in  applying  the  technique  included  variations  in  window  lengths,  both  in  the  time 
domain  over  which  the  ultrasonic  signal  was  selected  and  in  the  frequency  domain  over 
which  the  spectral  ratio  slope  was  taken.  A  conservative  estimate  of  error  in  applying  the 
technique  is  ±  20%. 


EXPERIMENTAL  RESULTS 
Hysteresis  Loop 

Table  1  contains  the  results  of  hysteresis  loop  attenuation  measurements  on  the  room 
dry  and  water-saturated  cylinder  of  Sierra  White  granite.  The  five  columns  list  the 
frequency,  strain  amplitude,  saturation,  Young's  modulus  in  GPa,  and  extensional 
attenuation  (1000/Qg).  The  strain  amplitude  is  the  maximum  peak  to  peak  amplitude  of  the 
recorded  displacement  signal  divided  by  the  ring  separation.  The  Young's  modulus  is 
calculated  from  the  corresponding  axial  stress  that  was  recorded. 

The  hysteresis  loop  for  the  dry  sample  is  shown  in  Figure  4.  In  this  plot  the  horizontal 
axis  is  strain  in  units  of  microstrain  and  the  vertical  axis  is  stress  in  units  of  MPa  times  ten. 
This  hysteresis  loop  was  averaged  from  128  cycles  of  sinusoidal  stress  applied  to  the 
sample  end  at  a  frequency  of  0.1  Hz.  The  time  sequence  of  stress  and  strain  follows  the 
hysteresis  loop  in  a  clockwise  rotation,  so  that  as  the  stress  is  reversed  the  strain  lags 
behind  and  forms  the  loop.  The  extensional  attenuation  factor  was  obtained  by  directly 
measuring  the  phase  angle  between  the  original  sinusoidal  signals  with  Fourier  analysis. 


Resonant  Bar 

Table  2  contains  the  result  of  resonant  bar  measurements  on  the  dry  and  saturated 
sample  of  Sierra  White  granite.  In  the  table  are  columns  of  center  resonant  frequency,  the 
harmonic  of  the  resonance,  the  saturation,  Young's  modulus  in  GPa,  and  extensional 
attenuation  (1000/Qg).  Young's  modulus  was  calculated  from  the  square  of  the  extensional 
velocity  times  the  density.  The  fundamental  resonance  peak  and  the  8th  harmonics  for  the 
dry  and  saturated  sample  are  shown  in  Fig.  5.  While  the  center  frequencies  of  the  two 
resonances  are  nearly  the  same,  the  much  greater  attenuation  for  the  saturated  sample  is 
indicated  by  the  broad  width  of  the  resonance.  The  higher  harmonics  tended  to  be 
distorted,  and  the  highest  frequency  for  the  water  saturated  measurements  was  limited  by 
this  distortion.  The  highest  frequency  for  the  dry  measurements  was  limited  by  the  upper 
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frequency  capacity  of  the  phase-lock  amplifier. 


Ultrasonic 

Table  3  contains  the  results  of  the  benchtop  ultrasonic  velocity  and  attenuation 
measurements  on  Sierra  White  granite  samples.  The  measurements  are  grouped  vertically 
in  the  table  according  to  the  spectral  bandwidth  of  the  ultrasonic  transducers,  with  the 
highest  frequency  (0.6  - 1.1  MHz)  at  the  top,  followed  by  the  intermediate  frequency  (100 
-  200  kHz),  and  lowest  frequency  (40  -  80  kHz)  at  the  bottom.  The  columns  identify  the 
spectral  bandwidth  of  the  paired  ultrasonic  transducers,  a  number  code  for  the  sample  that 
is  explained  at  the  table  bottom,  the  saturation,  P  and  S  wave  velocities  in  km/s,  Young's 
modulus  in  GPa,  and  the  compressional  (P),  shear  (S),  and  extensional  (e)  attenuations 
(1000/Q).  The  equivalent  ultrasonic  Qe"1  was  computed  in  order  to  compare  the  Qp1  and 
Qs1  measured  in  the  ultrasonic  wave  propagation  experiments  with  the  Qe'1  measured  in 
the  resonant  bar  and  hysteresis  loop  experiments.  If  the  material  is  assumed  to  be  isotropic 
and  the  attenuation  is  small,  then  complex  moduli  can  be  substituted  for  real  moduli  in 
linear  elasticity  relations  and  the  following  relationship  follows  (Winkler  and  Nur,  1979): 

(I-v)  (l-2v)  Qp'1  =  (1-n)  <^'*  -  2v  (2-v)  1 

where  v  is  Poisson's  ratio. 

P  and  S  wave  velocities  and  attenuations  were  also  measured  for  the  smaller  0.15  m 
Sierra  White  granite  sample  as  a  function  of  oo.. fining  pressure  for  both  the  dry  and  water 
saturated  conditions.  The  P  and  S  wave  velocity  data  are  presented  in  Figure  7  and  the 
calculated  P,  S,  and  e  attenuations  (1000/Q)  are  shown  in  Figure  8. 

Discussion 

The  experimental  data  on  Sierra  White  granite  determined  with  the  hysteresis  loop, 
resonant  bar,  and  ultrasonic  techniques  provide  a  good  basis  for  observations  and 
comparisons  concerning  the  effects  of  frequency,  saturation,  pressure,  and  experimental 
technique  on  attenuations  and  moduli.  Multiple  experiments  were  conducted  on  oriented 
samples  prepared  from  the  same  original  block  of  Sierra  White  granite.  A  broad  range  of 
frequencies  from  0. 1  Hz  to  1  MHz,  w '•rc  investigated,  and  there  was  partial  overlap  of 
frequencies  between  100  to  200  kHz  for  resonant  bar  and  ultrasonic  measurements.  Dry 
and  water  saturated  results  were  obtained  with  all  of  the  techniques.  Finally,  the  pressure 
dependence  of  attenuation  to  100  MPa  was  also  determined  with  dry  and  saturated 
ultrasonic  measurements  in  the  0.6  - 1.1  MHz  frequency  band. 
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The  attenuation  and  moduli  data  for  the  dry  samples  in  Tables  1, 2,  and  3  show 
surprisingly  good  agreement  for  Young's  modulus  but  a  complicated  result  for  attenuation. 
The  Young's  modulus  measurements  range  from  33.8  to  35.2  GPa  for  the  hysteresis  loop, 
the  fundamental  resonance  of  the  resonant  bar,  and  the  ultrasonic  P  and  S  wave 
measurements  over  two  frequency  bands  (100-200  kHz  and  0.6  - 1.1  MHz).  No 
monotonic  trend  exists  as  a  function  of  frequency,  and  the  4%  variation  can  be  considered 
as  part  of  experimental  error.  Although  desirable,  modulus  dispersion  over  the  broad 
frequency  range  cannot  be  accurately  established  with  this  data  due  to  limitations  of 
experimental  precision  and  sample  variability. 

The  non-ultrasonic  extensional  attenuations  are  comparable  for  the  dry  samples,  with 
the  hysteresis  loop  Qg*1  equal  to  0.0083  and  the  resonant  bar  value  equal  to  0.00526  at  the 
fundamental  resonance.  Between  the  fundamental  (10.436  kHz)  and  the  3rd  harmonic  of 
the  fundamental  (39.894  kHz)  resonant  bar  attenuation  monotonically  increases  to 
0.00662,  then  at  higher  harmonics  only  varies  between  0.00616  and  0.00795.  The 
corresponding  Young's  modulus  data  for  the  harmonics  decrease  between  the  fundamental 
and  3rd  harmonic,  and  thereafter  shows  random  variations  between  31.5  and  33.6  GPa. 
Some  of  the  negative  modulus  dispersion  between  the  fundamental  and  3rd  harmonic  may 
be  due  to  geometrical  effects,  which  would  also  contribute  to  higher  Qg'1,  but  it  is  not  clear 
why  the  trend  apparently  stops  at  the  3rd  harmonic. 

In  contrast  with  the  hysteresis  loop  and  resonant  bar  results,  the  ultrasonic  experiments 
show  that  there  is  a  strong  dependence  of  attenuation  on  frequency  (Table  3).  For  the  dry 
Sierra  White  granite  samples  the  three  ultrasonic  P-wave  attenuations  (1000/Qp)  at  0.6  - 1.1 
MHz,  100-200  kHz,  and  40-80  kHz  were  0.065, 0.028,  and  0.020,  respectively. 
Therefore  the  P-wave  attenuation  in  the  highest  frequency  band  was  over  three  times 
greater  than  in  the  lowest  frequency  band.  Similarly,  S-wave  attenuation  (1000/Qs)  in  the 
0.6  - 1.1  MHz  frequency  band  was  0.047,  or  approximately  four  and  one-half  times 
greater  than  the  0.01 1  value  measured  in  the  100-200  kHz  band.  The  calculated 
extensional  attenuations  (1000/Qe)  were  0.058  and  0.024  in  the  same  corresponding 
frequency  bands. 

The  frequency  dependence  of  ultrasonic  attenuation  raises  a  number  of  questions. 

First,  it  obviously  negates  the  assumption  in  the  spectral  ratio  calculation  that  attenuation  is 
constant  over  the  spectral  bandwidth  of  the  pulse.  Therefore  the  given  attenuations  must 
be  treated  as  estimates  for  the  limited  bandwidths  over  which  the  spectral  ratio  was  taken. 
Second,  the  trend  implies  that  frequency-dependent  scattering  is  probably  an  important 
attenuation  mechanism  for  ultrasonic  pulses  at  these  frequencies. 

A  particularly  interesting  feature  of  the  ultrasonic  and  resonant  bar  results  for  the  dry 
samples  is  the  overlap  in  frequency  from  approximately  100  to  200  kHz  (Tables  2  and  3). 
As  discussed  previously,  the  Young's  modulus  data  from  the  two  techniques  are  equivalent 
to  within  experimental  error.  The  calculated  ultrasonic  Qg'1  of  0.024  is,  however, 
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approximately  three  to  four  times  larger  than  the  resonant  bar  Qg*1,  which  range  from 
0.0065  to  0.008  for  the  higher  harmonics  in  the  same  frequency  band.  In  addition, 
resonant  bar  attenuations  to  the  20th  harmonic  (212  kHz)  on  the  dry  sample  show  no  clear 
trend  for  attenuation  changing  with  frequency.  Therefore,  a  fundamental  disagreement 
exists  between  the  attenuations  measured  with  the  ultrasonic  and  resonant  bar  techniques. 
The  difference  may  decrease  for  the  water  saturated  samples.  Although  ultrasonic 
measurements  on  saturated  samples  were  not  made  with  the  100-200  kHz  transducers,  the 
difference  would  probably  be  reduced  by  the  proportionally  larger  increase  in  resonant  bar 
attenuation  upon  saturation.  The  difference  may  also  decrease  at  lower  frequencies.  The 
trend  for  dry  P-wave  ultrasonic  attenuations  is  down  for  decreasing  frequencies  (Table  3), 
and  at  some  frequency  lower  than  50  kHz  the  resonant  bar  and  pulse  attenuations  may 
agree.  The  reasons  for  the  difference  between  ultrasonic  and  resonant  bar  attenuations  are 
not  entirely  clear.  One  possibility  is  that  the  scattering  attenuation  mechanism  that  elevates 
ultrasonic  pulse  attenuations  does  not  operate  equally  on  the  standing  waves  in  an 
extensional  resonant  bar  experiment. 

The  effect  of  water  saturation  on  attenuation  and  moduli  in  Sierra  White  granite  also 
contrasts  the  hysteresis  loop  and  resonant  bar  results  with  the  ultrasonic  pulse  results.  In 
the  hysteresis  loop  experiment,  water  saturation  caused  Young's  modulus  to  decrease  from 
33.8  to  31.2  GPa  and  Qg"1  to  increase  dramatically  from  0.0083  to  0.026  (Table  1). 
Similarly,  in  the  resonant  bar  experiment,  water  saturation  caused  a  large  increase  of  Qg'1 
from  0.00526  to  0.0397  at  the  fundamental  resonance  and  only  a  slight  increase  in  Young’s 
modulus,  from  35.2  to  35.3  GPa,  a  change  that  is  within  the  experimental  error.  Both 
attenuation  and  Young's  modulus  decreased  at  higher  harmonic  resonances,  however,  with 
attenuation  dropping  to  0.0255  and  Young's  modulus  to  31.8  GPa  at  the  8th  harmonic.  In 
summary,  saturation  causes  a  slight  decrease  in  Young's  modulus  at  the  lower  frequency, 
no  change  at  the  fundamental  (10.4  kHz)  frequency,  and  increases  by  a  factor  of  three  and 
a  factor  of  seven  and  one-half  the  extensional  attenuation  for  the  hysteresis  loop  and 
resonant  bar  attenuations,  respectively. 

The  large  attenuations  associated  with  saturation  can  be  interpreted  as  a  results  of 
viscous  fluid  flow  within  thin  microcracks  and  pores  of  the  granite.  Two  dissipation 
mechanisms  are  possible  with  fluid  flow.  The  rest  is  inertial  shear  between  crack  surfaces 
and  bulk  pore  fluid  (Biot,  1 956a, b).  The  energy  loss  associated  with  this  mechanism, 
however,  is  very  small  in  a  low-porosity  granite  at  sub-ultrasonic  frequencies.  The  second 
mechanism  is  viscous  loss  within  pore  fluid  due  to  pressure  gradients  between  adjacent 
pore  volumes  induced  by  dynamic  elastic  deformation.  This  mechanism  has  received 
considerable  attention  (O'Connell  and  Budiansky,  1977;  Mavko  and  Nur,  1979;  Murphy, 
1982;  Bulau  et  al.,  1984)  but  specific  models  of  pore  distributions  and  pressure  gradients 
have  not  been  developed. 
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The  ultrasonic  pulse  results,  on  the  other  hand,  probably  combine  fluid  flow  losses 
with  scattering.  In  the  0.6  -  1.1  MHz  frequency  band  Qs_1  and  Qg"1  increase  upon 
saturation  (Table  3),  although  not  nearly  as  much  as  at  lower  frequencies.  This  can  be 
interpreted  as  a  fluid  flow  related  attenuation  mechanism.  The  Qp"1,  however,  actually 
decreases  upon  saturation,  an  observation  that  correlates  well  with  a  scattering  mechanism. 
Winkler  (1983)  has  concluded  that  scattering  is  quite  prevalent  in  sandstones  at  ultrasonic 
frequencies.  His  calculations  of  the  scattering  cross-section  of  a  dry  and  saturated  single 
spherical  scatterer  indicate  that  the  dry  case  will  results  in  approximately  twice  the  scattering 
loss  of  a  fluid-filled  scatterer.  The  measured  reduction  of  Qp"1  upon  saturation  agrees  with 
the  sign  of  the  calculation  but  not  quite  the  magnitude,  probably  because  scattering  is 
combined  with  viscous  fluid  loss.  In  the  40-80  kHz  frequency  band  the  reverse  effect  is 
observed.  The  water  saturated  Qp"1  (0.064)  is  much  larger  than  the  dry  (0.020).  At  these 
lower  frequencies  fluid  flow  loss  may  dominate  over  scattering  loss. 

Confining  pressure  closes  or  reduces  the  size  of  microcracks,  increasing  ultrasonic 
velocities  and  reducing  attenuations  (Figures  7  and  8).  The  dry  and  saturated  properties 
also  converge  with  increasing  pressure.  The  microcrack  population  responsible  for  the  P 
wave  scattering  is  apparently  removed  by  the  application  of  45  MPa  of  effective  confining 
pressure.  The  convergence  of  dry  and  saturated  attenuations  also  indicates  that  the  fluid 
flow  mechanism  is  reduced  as  thin  microcracks  close  under  pressure.  At  the  highest 
confining  pressure  the  minimum  extensional  attenuation  is  0.022,  which  is  still  a  factor  of  3 
to  4  times  greater  than  the  dry  hysteresis  loop  and  resonant  bar  attenuations.  Ultrasonic 
scattering  may  therefore  still  occur  at  high  pressures  because  microcracks  do  not  close 
entirely. 

The  experimental  data  in  this  study  is  graphically  summarized  in  Figures  9  and  10  along 
with  several  sets  of  previous  experimental  results.  These  plots  give  the  extensional  mode 
attenuations  (1000/Qg)  versus  log  frequency  for  dry  (Fig.  9)  and  saturated  (Fig.  10)  Sierra 
White  granite  obtained  in  this  study  plus  pertinent  results  on  various  other  granites.  The 
hysteresis  loop,  resonant  bar,  and  ultrasonic  pulse  data  are  labelled  in  each  plot  The 
results  of  others  are  number-coded  according  to  the  key  in  the  figure  caption.  Spencer 
(1981)  and  Liu  and  Peselnick  (1983)  conducted  experiments  with  the  hysteresis  loop 
technique  on  Oklahoma  and  Westerly  granites,  respectively.  Spencer  collected  data  on  a 
vacuum  dry  sample  as  a  function  of  frequency  from  approximately  5  Hz  to  400  Hz,  and  the 
solid  line  in  Fig.  9  represents  an  upper  bound  for  the  data  he  collected.  The  Qg"1  reported 
by  Liu  and  Peselnick  (1983)  on  a  room  dry  sample  have  been  increased  by  a  factor  of  4 
(discussed  in  the  first  section  of  this  report).  Also  included  in  these  plots  are  resonant  bar 
results  by  Murphy  (1984)  and  Winkler  et  al.  (1979)  on  Sierra  White  granite,  Gordon  and 
Davis  (1968)  on  Tupper  Lake  granite,  Johnston  and  Toksoz  (1980)  and  Tittmann  (1984) 
on  Westerly  granite.  These  attenuation  data  are  all  Qg"1  except  for  Tittmann's  results, 
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which  are  Qs_1  collected  in  the  flexural  and  torsional  modes. 

The  attenuation  data  in  Fig.  9  for  room  dry  granite  specimens  are  fairly  well  clustered 
between  0.005  and  0.010  (Qg*1)  at  frequencies  from  0.1  Hz  to  the  highest  frequency 
resonant  bar  experiment  (212  kHz).  These  data  include  the  hysteresis  loop  and  resonant 
bar  experiments,  plus  the  results  of  Liu  and  Peselnick  (1983),  Tittmann  (1984), 

Winkler  et  al.  (1979),  and  Johnston  and  Toksoz  (1980)  on  Westerly  and  Sierra  White 
granites.  Collectively,  the  data  indicate  a  lack  of  frequency  dependence  for  attenuation  in 
these  rom  dry  granites,  and  surprisingly  good  agreement  between  the  hysteresis  loop  and 
resonant  bar  measurements. 

Attenuation  in  nominally  "dry"  rock  is  closely  related  to  moisture  content  and 
microcrack  population.  An  experimentally  measured  Q*1  of  0.1  for  a  basalt  exposed  to 
high-humidity  air  can  be  reduced  to  2*10*4  by  following  a  strong  outgassing  procedure 
Tittmann  (1977).  In  addition,  Spencer  (1981)  showed  that  hysteresis  loop  attenuation 
measurements  for  vacuum  dry  samples  are  similarly  very  low.  The  line  segment  associated 
with  Spencer  (1981)  in  Fig.  9  is  actually  an  upper  bound  for  many  measurements  on 
Oklahoma  granite  that  were  much  closer  to  zero  attenuation.  The  severe  attenuation  caused 
by  small  amounts  of  surface  adsorbed  water  is  probably  a  combination  of  viscoelastic  and 
electrostatic  responses  to  stress  by  thin  water  layers  in  the  microcracks  of  what  are 
otherwise  "dry"  rocks  (Bulau  et  al.,  1984).  Therefore,  the  clustering  of  Qg'1  factors 
between  0.005  and  0.010  is  a  reflection  of  typical  microcrack  populations  for  these  granites 
and  laboratory  ambient  humidities. 

The  attenuation  factors  for  the  water-saturated  granites  (Fig.  10)  are  substantially  larger 
than  for  the  dry  samples  and  exhibit  a  complicated  dependence  on  frequency.  A  Qg*1  of 
0.026  was  measured  by  the  hysteresis  loop  experiments  at  0.1  Hz,  approximately  three 
times  greater  than  for  the  dry  sample.  A  greater  increase  was  also  measured  and  a  resonant 
bar  experiments,  although  between  10  kHz  and  90  kHz  Qg’1  is  decreasing  as  a  function  of 
increasing  frequency. 

The  experiments  on  Sierra  White  granite  reported  in  this  paper  are  consistent  with,  but 
do  not  fully  confirm,  a  broad  relaxation  peak  in  the  extensional  mode  attenuation  spectrum 
for  saturated  rocks  reported  by  Spencer  (1981).  The  dashed  line  in  Fig.  10  outlines  the 
shape  of  the  relaxation  peak  with  a  center  frequency  of  around  1  kHz.  The  flanks  of  the 
peak  are  partially  defined  by  Spencer's  (1981)  hysteresis  loop  data  between  5  and  400  Hz 
for  Oklahoma  granite,  and  the  resonant  bar  data  between  10  and  90  kHz  for  Sierra  White 
granite.  Murphy's  (1984)  resonant  bar  data  on  Sierra  White  granite  apparently  fall  close 
enough  to  the  relaxation  peak  so  as  to  not  exhibit  dependence  of  attenuation  on  frequency. 
The  hysteresis  loop  data  at  0.1  Hz  is  also  consistent  with  the  relaxation  peak,  but  the 
ultrasonic  Qg*1  may  be  high  due  to  scattering.  The  Davis  and  Gordon  (1968)  data  point  on 
saturated  Tupper  Lake  granite  is  somewhat  high,  as  it  was  for  the  dry  sample,  but  still 
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roughly  follows  a  trend  of  decreasing  attenuation  at  the  higher  frequencies.  The  only 
serious  discrepancy  is  found  in  the  torsional  and  flexural  attenuations  on  Westerly  granite 
by  Tittmann  (1984),  which  are  well  below  all  of  the  other  data.  The  reason  for  this  is  not 
clear,  although  shear  attenuation  should  be  similar  to  or  greater  than  extensional  mode 
attenuation  (Murphy,  1984;  Winkler  and  Nur,  1979). 

The  relaxation  peak  and  elevated  attenuation  for  the  saturated  granites  can  be  generally 
explained  by  attenuation  mechanisms  based  on  viscous  fluid  flow  within  the  microcracks  in 
response  to  a  propagating  stress  wave  (O'Connell  and  Budiansky,  1977;  Mavko  and  Nur, 
1979).  An  additional  attenuation  mechanism  suggested  by  Spencer  (1981),  based  on  his 
analysis  of  the  effects  of  various  saturants  on  the  relaxation  peak  and  Young's  modulus,  is 
the  dynamic  response  of  chemical  bonds  between  saturating  fluid  and  mineral  surfaces. 
Either  dissipation  process  results  in  a  peak  in  the  attenuation  spectrum,  with  a  transition 
from  lower  moduli  at  frequencies  below  the  center  frequency  to  higher  unrelaxed  moduli  at 
frequencies  above  the  center  frequency.  The  resonant  bar  results  in  this  study,  however, 
show  that  saturated  Young's  moduli  are  not  much  larger  than  the  saturated  hysteresis  loop 
value,  and  not  nearly  that  predicted  by  Spencer  (1981)  in  this  Cole-Cole  distribution  of 
relaxation  times.  One  possible  explanation  is  that  the  unjacketed  resonant  bar  may  have 
only  been  partially  saturated  during  the  saturated  experiment,  resulting  in  a  lower  Young's 
modulus. 

Further  experimental  investigations  utilizing  all  three  techniques  are  warranted  in  order 
to  fully  document  and  interpret  attenuation  in  granite.  In  this  study  the  ultrasonic  pulse 
technique  was  applied  in  three  frequency  bands,  from  100  to  200  kHz  and  0.6  to  1.1  MHz 
for  both  P  and  S  waves,  and  from  40  to  80  kHz  for  the  P  wave.  Although  ultrasonic 
attenuations  are  large,  particularly  for  dry  samples,  attenuation  is  clearly  decreasing  at 
lower  frequencies.  A  suitable  S  wave  transducer  in  the  40-80  kHz  band  needs  to  be 
developed,  and  dry  and  water  saturated  data  at  these  lower  frequencies  are  desirable.  At 
some  frequency  the  unrelaxed  ultrasonic  modulus  and  attenuation  should  decrease  and 
agree  with  the  hysteresis  loop  and  resonant  bar  data.  An  understanding  of  the  scattering 
attenuation  mechanism  and  attenuation  spectra  for  these  ultrasonic  pulses  in  intact  rock  may 
provide  some  insight  into  the  larger-scale  scattering  phenomena  that  occur  at  seismic 
frequencies  in  the  earth's  crust. 

The  relaxation  peak  that  is  postulated  to  occur  around  1  kHz  for  the  saturated  granite 
requires  further  definition.  The  results  of  Spencer  (1981)  and  Bulau  et  al.  (1984)  show  an 
increase  of  attenuation  towards  1  kHz;  the  data  collected  on  Sierra  White  granite  displays  a 
corresponding  decrease  in  the  resonant  bar  data  at  frequencies  from  10  to  90  kHz.  The 
resonant  bar  technique  may  be  used  for  a  complete  frequency  sweep  through  the  relaxation 
peak  by  reducing  the  fundamental  resonance,  either  through  loading,  in  the  manner  of 
Tittmann  (1984),  or  with  longer  samples.  In  order  to  obtain  fundamental  resonant 
frequency  of  around  1  kHz  samples  on  the  order  of  1.7  meters  long  would  have  to  be 
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prepared.  At  low  frequencies,  on  the  opposite  side  of  the  relaxation  peak,  the  hysteresis 
loop  technique  can  be  performed  over  a  broader  frequency  band,  in  the  manner  developed 
by  Spencer  (1981).  This  would  fill  in  the  frequency  dependence  over  the  seismic 
frequency  band  and  better  define  the  relaxation  peak. 


Conclusions 

In  this  study  a  comparison  has  been  made  of  the  moduli  and  attenuations  for  dry  and 
saturated  Sierra  White  granite  determined  with  the  hysteresis  loop,  resonant  bar,  and 
ultrasonic  techniques.  The  hysteresis  loop  and  resonant  bar  results  indicate  that  for  room- 
dry  granite  the  Young's  modulus  and  extensional  attenuation  are  nearly  constant  over  a 
broad  frequency  range  from  0.1  Hz  to  at  least  200  kHz.  Water  saturation  caused  a  dramatic 
increase  in  extensional  attenuation,  by  factors  of  four  to  seven  and  one-half,  but  had  a 
minimal  effect  on  Young's  modulus.  In  addition,  resonant  bar  attenuations  in  the  range 
10-90  kHz  decreased  with  increasing  frequency,  consistent  with  a  relaxation  peak  in  the 
extensional  attenuation  spectra  at  a  frequency  below  10  kHz.  The  attenuation  mechanism  is 
viscous  fluid  flow  within  the  microcracks  of  the  granite. 

Ultrasonic  attenuations  and  moduli  in  three  different  frequency  bands  show  that 
scattering  is  probably  an  important  attenuation  mechanism  at  these  frequencies.  Although 
the  dry  Young's  moduli  derived  from  P  and  S  wave  velocities  agree  with  the  resonant  bar 
and  hysteresis  loop  moduli,  in  the  100-200  kHz  frequency  band  the  ultrasonic  Qg*1  is 
larger  than  the  resonant  bar  values  by  a  factor  of  3  to  4.  Scattering  is  also  manifested  in  the 
effects  of  saturation  on  attenuation.  The  water  saturated  Qp*1  at  0.6  - 1.1  MHz  is  less  than 
dry  because  the  effective  cross-sectional  area  of  the  scattering  cracks  is  probably  reduced 
by  the  presence  of  fluid.  The  scattering  mechanism  also  combines  with  fluid  losses, 
however,  because  Qg'1  and  Qs-1  increase  in  the  saturated  sample  at  0.6  - 1.1  MHz. 
Confining  pressure  closes  or  reduces  the  size  of  microcracks,  increasing  velocities  and 
reducing  attenuation.  At  high  pressures  the  dry  and  saturated  velocities  and  attenuations 
converge. 

The  three  experimental  techniques,  hysteresis  loop,  resonant  bar,  and  ultrasonic,  each 
have  characteristics  that  can  be  uniquely  applied  to  various  requirements  for  measuring  and 
documenting  dynamic  properties  of  rock.  From  the  seismological  perspective  the  most 
desirable  and  relevant  laboratory  determinations  would  utilize  hysteresis  loops  because  the 
frequency  bandwidths  are  comparable.  Since  the  cyclic  stress  can  be  varied  to  any  desired 
amplitude,  large  strain-amplitude  measurements  can  also  be  made,  even  to  near  failure.  In 
addition,  large  samples  can  be  used  for  investigating  larger  scale  inhomogeneities  such  as 
fractures.  On  the  other  hand,  the  technique  is  experimentally  difficult  and  loses  resolution 
at  very  low  attenuations.  Ultrasonic  attenuations  are  probably  dominated  by  scattering 
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mechanisms.  Although  the  exact  attenuation  values  may  not  be  particularly  relevant  to 
lower  frequencies,  scattering  phenomena  at  lower  frequencies  may  perhaps  be  studied. 
Resonant  bar  experiments  have  the  greatest  resolution  at  very  low  attenuations.  In 
addition,  they  may  be  conducted  under  confining  pressure  and  with  different  saturations. 

In  addition  to  the  extensional  resonance  experiments  performed  in  this  study, 
complimentary  torsional  or  flexural  shear  measurements  may  also  be  made.  Resonant  bar 
measurements  are  probably  a  suitable  compromise  for  dry  measurements,  given  the  lack  of 
frequency  dependence  for  dry  samples  below  ultrasonic  frequencies. 
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TABLE  1 .  Attenuation  measurements  in  Sierra  White  granite  utilizing  the  hysteresis 
loop  technique. 


FREQUENCY 

STRAIN  SATURATION 

AMPLITUDE 

MICROSTRAIN 

YOUNG’S 

MODULUS 

(GPa) 

ATTENUATION 

(lOOCVOg) 

0.1  Hz 

0.50 

Dry 

33.8 

8.3 

1.92 

Water 

31.2 

26.0 
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TABLE  2.  Attenuation  measurements  in  Sierra  White  granite  utilizing 
the  resonant  bar  technique. 


FREQUENCY 

(kHz) 

HARMONIC 

SATURATION 

YOUNGS 

MODULUS 

(GPa) 

ATTENUATION 

(1000/Qe) 

10.436 

0 

Dry 

35.2 

5.26 

20.600 

1 

Dry 

34.3 

5.65 

30.462 

2 

Dry 

33.3 

6.11 

39.894 

3 

Dry 

32.2 

6.62 

50.064 

4 

Dry 

32.4 

6.16 

59.843 

5 

Dry 

32.2 

6.29 

69.135 

6 

Dry 

31.5 

6.70 

80.048 

7 

Dry 

32.4 

6.56 

90.245 

8 

Dry 

32.5 

6.50 

101.252 

9 

Dry 

33.2 

7.15 

110.824 

10 

Dry 

32.8 

6.70 

121.804 

11 

Dry 

33.3 

7.38 

132.567 

12 

Dry 

33.6 

6.61 

141.595 

13 

Dry 

33.1 

7.45 

152.075 

14 

Dry 

33.2 

7.95 

161.999 

15 

Dry 

33.2 

6.80 

171.243 

16 

Dry 

32.8 

7.54 

181.654 

17 

Dry 

32.9 

7.31 

192.130 

18 

Dry 

33.1 

7.00 

102.797 

19 

Dry 

32.9 

6.72 

211.928 

20 

Dry 

32.9 

6.51 

10.455 

0 

Water 

35.3 

39.7 

20.534 

1 

Water 

34.1 

35.7 

30.435 

2 

Water 

33.3 

35.9 

40.902 

3 

Water 

33.8 

29.1 

49.892 

4 

Water 

32.2 

30.7 

59.251 

5 

Water 

31.5 

26.8 

68.478 

6 

Water 

30.9 

26.8 

70.050 

7 

Water 

31.6 

23.0 

89.262 

8 

Water 

31.8 

25.5 
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TABLE  3.  Attenuation  measurements  in  Sierra  White  granite  utilizing  the  ultrasonic 
pulse  technique. 


FREQ. 

SAM. 

SAT. 

VELOCITY 

(km/sec) 

YOUNGS 

MODULUS 

(GPa) 

ATTENUATION 

(1000/Q) 

P  WAVE 

S  WAVE 

P 

S 

E 

0.6- 1.1 

1 

Dry 

3.686 

2.434 

34.8 

65 

47 

58 

MHz 

Water 

5.444 

3.127 

64.7 

46 

87 

75 

100-200 

2 

Dry 

3.593 

2.454 

33.8 

28 

11 

24 

kHz 

40-80 

3 

Dry 

3.536 

20 

kHz 

Water 

5.331 

64 

Samples:  1  =  0. 152  m  long  cylinder,  2  =  0.203  m  thick  block,  3  =  0.508  m  long  cylinder 
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Peizoelectric  Crystal 


Tuff  Specimen 


Peizoelectric  Crystal 


Figure  2.  Schematic  diagram  of  the  benchtop  resonant  bar  apparatus. 
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Figure  3.  Schematic  diagram  of  the  experimental  apparatus  for  ultrasonic  velocity  and  attenuation 
measurements  under  pressure. 
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SIERRA  WHITE  GRANITE 


0  0.2  0.4  0.6 

AXIAL  MICROSTRAIN 


Figure  4.  Hysteresis  loop  for  dry  Sierra  White  granite  at  a  strain  amplitude  of  0.56 
microstrain.  The  vertical  scale  is  axial  stress  in  MPa.  The  attenuation  factor 
calculated  for  this  sample  is  0.0083. 
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Figure  5.  Fundamental  (left)  and  8th  harmonic  (right)  extensional  resonance  for  room-dry  (top) 
and  water  saturated  (bottom)  Sierra  White  granite. 
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Figure  6.  Ultrasonic  (0.6-1. 1  MHz)  P-  and  S-wavc  signals  propagated  through  the  0.152  m 
sample  of  dry  (upper  left)  and  water  saturated  (upper  right)  Sierra  White  granite  and  through  the 
0.15  m  aluminum  standard  (bottom). 
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Figure  7.  P  and  S  wave  velocities  for  dry  and  water  saturated  Sierra  White  granite  measured  as  a 
function  of  confining  pressure. 
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Figure  8.  Ultrasonic  (0.6  - 1.1  MHz)  P  and  S  wave  attenuation,  and  calculated  extensional 
attenuation  (1000/Q),  for  dry  and  water  saturated  Sierra  White  granite  measured  as  a  function  of 
confining  pressure. 
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GRANITES 


ATTENUATION  VERSUS  LOG  FREQUENCY 


LOG  FREQUENCY,  Hz 


Figure  9.  Experimental  attenuation  measurements  versus  log  frequency  for  dry  granite. 
Hysteresis  loop,  resonant  bar,  and  ultrasonic  pulse  results  are  labelled.  Other  results: 

(1)  Liu  and  Peselnick,  1983;  (2)  Spencer,  1981;  (3)  Winkler  et  al.,  1979;  (4)  Tittmann,  1984; 
(5)  Johnston  and  Toksoz,  1980;  (6)  Gordon  and  Davis,  1968 
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ATTENUATION  (1000/Q) 


GRANITES 


ATTENUATION  VERSUS  LOG  FREQUENCY 


LOG  FREQUENCY,  Hr 


Figure  10.  Experimental  attenuation  measurements  versus  log  frequency  for  water  saturated 
granites.  Hysteresis  loop,  resonant  bar,  and  ultrasonic  pulse  results  are  labelled.  Other  results: 
(1)  Spencer,  1981;  (2)  Murphy,  1984;  (3)  Gordon  and  Davis,  1968;  (4)  Tittmann,  1984; 
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